Evolving fracture patterns: columnar joints, mud cracks, and polygonal terrain 
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When cracks form in a thin contracting layer, they sequentially break the layer into smaller and smaller 
pieces. A rectilinear crack pattern encodes information about the order of crack formation, as later cracks tend 
to intersect with earlier cracks at right angles. In a hexagonal pattern, in contrast, the angles between all cracks 
at a vertex are near 120°. However, hexagonal crack patterns are typically only seen when a crack network 
opens and heals repeatedly, in a thin layer, or advances by many intermittent steps into a thick layer. Here it 
is shown how both types of pattern can arise from identical forces, and how a rectilinear crack pattern evolves 
towards a hexagonal one. Such an evolution is expected when cracks undergo many opening cycles, where 
the cracks in any cycle are guided by the positions of cracks in the previous cycle, but when they can slightly 
vary their position, and order of opening. The general features of this evolution are outlined, and compared 
to a review of the specific patterns of contraction cracks in dried mud, polygonal terrain, columnar joints, and 
eroding gypsum-sand cements. 



I. INTRODUCTION 



Cracks in cooling or drying media can form captivating patterns of connected networks, such as the 
artistic craqulere patterns sometimes seen in pottery glazes, to those found in dried mud, or the polygonal 
networks covering the polar regions of Earth and Mars. Two types of pattern are common. The first is a 
rectilinear pattern, such as may form as an homogeneous slurry is dried or cooled uniformly ITl l2T1l6T1l . 
The verticies of these patterns typically contain one crack path intersecting another straight crack at right 
angles, a so-called T-junction. Secondly are crack networks with a more regular hexagonal pattern, such as 
the famous columnar joints of the Giant's Causeway (8j [3T] . These are dominated by cracks intersecting 
at 120°, or Y-junctions. Curiously, examples of both extreme types of patterns can sometimes be found in 
the same systems, such as garden-variety mud-cracks. These, and other examples, are shown in Fig. [T] In 
general, hexagonally ordered patterns may form as the result of either desiccation or thermal contraction 
cracks, and there is a deep physical symmetry between these two mechanisms (5JI53). Furthermore, they 
may be roughly classified into crack patterns that have evolved in both time and space, such as a network 
of crack tips that delimit the prismatic forms of columnar joints as they advance l43ll . and those that have 
recurred in the same place, but evolved over time as the pattern repeatedly healed, and re-cracked l2T1l35ll . 

Here it is shown that both hexagonal and rectilinear crack patterns can arise naturally from the same 
simple assumptions of fracture behaviour, and how a rectilinear, T-junction dominated pattern can develop 
into to a hexagonal pattern, with Y-junctions. This paper is roughly divided into two parts. In the first part 
(Sections 2,3), the relevant physics of linear elastic fracture mechanics in a thin brittle layer is summarised, 
and used to develop a general model for the evolution of a rectilinear crack pattern into a hexagonal one, 
through the cyclic opening of cracks. The second part (Sections 4-7) is a review of the ordering of crack 
patterns in a number of geophysical systems: desiccation cracks in clays that are repeatedly wetted and 
dried; thermal contraction cracks in permafrost, also known as polygonal terrain; columnar joints in lava 
and starch; and cracks in eroding gypsum sand layers. These patterns will each be explored, in turn, and 
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FIG. 1: Hexagonal crack patterns are found naturally in many settings. They can be driven by thermal or desiccation 
contraction, and appear to require the opportunity to evolve, either in space or time. Shown here are examples of colum- 
nar joints in (a) lava at the Giant's Causeway and (b) in desiccated starch, (c) desiccation cracks on the cemented sands 
on the windward slopes of dunes (reproduced from Ref. 1 10|), (d) 'young' thermal contraction cracks, or polygonal 
terrain, in Taylor Valley, Antarctica, which is speculated to be in the process of evolving into a mature pattern such 
as (e) that in Beacon Valley, Antarctica, and desiccation-induced mud-crack patterns (f) in the lab, and (g) as found 
naturally in soil. Images (a,d,e) were taken obliquely, and the scale estimated by feature sizes. In (b) the starch columns 
are shown inverted, with the original drying surface indicated by a dashed line 



shown to obey the same general physical mechanism of crack pattern evolution. In the case of desiccation 
cracks in clays new experiments are also reported, which show that the length-scale of the evolution of 
a rectilinear mud-crack pattern towards a hexagonal one scale linearly with the thickness of the cracking 
layer. 



II. CONTRACTION CRACKS IN THIN LAYERS 



When an elastic body cools or dries it tends to shrink, and may crack. Here the case of a brittle film 
under tensile stress, and firmly attached to a much thicker, non-cracking substrate, as sketched in Fig. [2ja), 
is reviewed. As a drying body looses water, it develops an internal pore pressure as a result of capillary 
effects. In the absence of cracks the equations of stress equilibrium and poro-elasticity are easily solved 
to give an in-plane stress, <xo, proportional to this pressure, in a drying film. Similarly, if such a film 
is cooled, it will have an in-plane stress proportional to the temperature change, and to the mismatch in 
thermal expansivities of the film and substrate. The similarity between the patterns of cracking during 
drying and cooling reflects the underlying similarity between the physics of poro-elasticity, and thermo- 



FIG. 2: Schematic of thin film fracture, (a) A film of height h under stress cr and with elastic properties E\ and V\ is 
adhered to a much thicker substrate with potentially different elastic properties E 2 and v 2 . (b) In such a film a channel 
crack may open, but not penetrate into the substrate, and grow across the film, releasing stress preferentially in the 
direction perpendicular to it's growth direction, (c) If many such cracks form, the cracks will interact with each other. 
Later cracks will curve to intersect earlier ones at right angles, the direction that maximises the strain energy release 
rate as they grow. 

elasticity l5ll53l. 

Tension is released when a crack opens in the brittle layer, as the normal stresses in the film must vanish 
along the new crack surfaces. The substrate will resist deformation, however, and exert a restoring force 
on the film- substrate interface, approximately proportional to the local in-plane strain (66] [68]|. For a long 
straight crack in a thin film, as shown in Fig. [2^b), these conditions give rise to an exponential dependence 
of the height-averaged film stress on position x away from the crack, 



where cr ± and cry are the stress normal and parallel to the crack face, and v\ is Poisson's ratio of the film 
1581 . The length-scale over which stresses relax, /, is proportional to the film thickness, h, 



The term g here is given numerically in reference 1 32 ] and depends on the Dundurs parameters a and ft, 
ratios involving the plane-strain elastic moduli E\ and E 2 and Poisson ratios v\ and v 2 of the film and 
substrate. If there is no elastic mismatch between the cracking layer and its substrate, / = 2.0 h, while for a 
compliant film on an un-yielding substrate / is between 1.2 and 1.3 h, depending on v\. In the limit of very 
soft (or fluid) substrates, / diverges, and the film behaves as an unconstrained layer. 

Once a crack has formed, the Griffith condition for fracture propagation compares the energy required 
to create new fracture surface, G c , to the change in elastic energy in the cracking body for an infinitesimal 
extension of the crack, G, and predicts crack growth when more energy would be released, then consumed 
l36l . While the first term, the critical energy release rate, is a material parameter, the second term depends 
on the geometry of the crack and its surroundings. For an isolated straight film crack that is long compared 
to the film thickness [ 32[[88), the strain energy release rate 



cr ± = (T (l - e~ x/ \ <r,| = <r (l - vie~ x/l ) 



(1) 



/ = -hg(a,0). 



(2) 




(3) 



Thus, when the film stress reaches some critical value cr c = y/EG c /l, the film may crack. Since G c is 
a material constant, thinner layers will require higher stresses to crack than thicker layers, a fact well- 
established in the engineering literature fTTll30ll83l . 
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As stress develops, it can be relieved by the sequential fragmentation of the film. This can be illustrated 
by considering an array of parallel cracks, each separated by a distance A that is large compared to /. If each 
crack relaxes stress as per Eqn. [T] then the energy release rate, per crack, of such an array is 

G = G 4 tanh (^)' (4) 

where cr c is the stress required to open a single isolated crack, or critical cracking stress l66l . As drying 
or cooling proceeds, the background stress cr increases, and the allowed crack spacing X will shrink. New 
cracks will thus be able to open, intermittently, breaking the layer into smaller and smaller pieces. Counter- 
intuitively, this process does not continue indefinitely, no matter how high the film stress gets. When the 
cracks are closer than a few film thicknesses to each other then stresses can no longer be approximated by a 
through-thickness-average stress. More accurate simulations of the stress state between two parallel cracks 
have shown that, for cracks closer than some critical separation A c ~ h, the surface stress in the region 
between the two cracks is compressive, and will tend to close any further flaws or potential cracks l3l[67l. 
For cracks in a compliant layer on a hard substrate, experiments have confirmed that the crack spacing is 
limited to about 3 times the film thickness (T]|2T][29l[54j[6Tl. It is likely that in other cases the saturation 
spacing scales with the stress relaxation length /, although this has not been rigorously demonstrated. If so, 
then where there is no difference in the elastic properties of the cracking layer, and its substrate, the limiting 
crack spacing should be approximately 5 times the depth of the cracks. 

When multiple cracks form they will interact with each other and pattern the film. Rather than the simple 
parallel cracks just considered, cracks will nucleate randomly on material defects or along the faces of pre- 
existing cracks and grow until they hit something, such as a boundary or other crack. Each crack will tend 
to grow in the direction that maximises G - G c , the difference between the strain energy released, and the 
fracture energy consumed (361 . This is a local principle, as the energy release rates are calculated at the 
point and time of the crack growth. Since, as shown in Eqn. [TJ any pre-existing crack is more efficient at 
relieving stress perpendicular to its growth direction, than parallel to it, a growing crack will tend to curve to 
intersect at right angles to it. As shown in Fig. |2jc), the resulting T-junctions thus encode time information 
that can be used to partially reconstruct the order in which cracks form Q. 

III. A MODEL FOR EVOLVING CRACK PATTERNS 

Cracks in a homogeneous film will tend to pattern the surface into roughly equally- sized rectilinear 
pieces, bounded with right-angled T-junctions. However, many natural crack patterns, such as columnar 
joints, show hexagonal planforms, and equi-angled Y-junctions. In both cases Euler's formula implies that 
the average number of neighbours for any polygon is 6, as long as cracks do not cross each other, or vertices 
form from more than three cracks. It is possible to continuously deform a rectilinear crack network into a 
hexagonal one (allowing vertices to pass each other, if necessary ED), if the pattern is allowed to evolve. 
Such an evolution can be explained as the result of meeting three conditions: (1) that cracks recur or advance 
repeatedly, (2) that the previous positions of cracks act as lines of weakness, guiding the next iteration of 
cracking; and (3) that the order of opening of cracks may change in each iteration. 

In order to develop a rectilinear network of cracks into a hexagonal one, the pattern must have some 
means of changing with time. Polygonal terrain consists of networks of thermal contraction cracks that 
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FIG. 3: Schematic model for the evolution of a rectilinear crack pattern into a hexagonal pattern, (a) In the first cycle, 
the order of crack opening is captured by the shape of the T-junction, (b) If the cracks heal and recur a new crack can be 
guided by the previous cycle, but deflected by the asymmetric stress concentration it feels near a vertex, in which case 
(c) the vertex is shifted, and the crack paths curved, (d) This process may repeat, and evolve the crack pattern. Over 
many cycles (e), with cracks equally likely to approach a vertex from any branch, symmetry will lead to a Y-junction. 
(f) Close to the vertex, however, the order of cracking in any particular cycle may still curve cracks slightly, (g) If the 
evolution of the crack network is stacked in space, it suggests a model for extended patterns such as columnar joints. 

open during the cold polar winters, and which thaw, and heal, during summers. The crack pattern forms 
anew each year. For columnar joints the pattern of growing cracks at any particular moment is confined to 
a thin region, approximating a 2D film. Each crack advances step-wise, and the prismatic columns are the 
forms delimited by the planar crack network as it advances through space. 

Even if a layer cracks and heals repeatedly, the pattern may not evolve. If the system has no memory 
of previous patterns, each cycle will lead to a novel, but statistically similar, rectilinear crack pattern. If 
there is memory of vertex positions only, perhaps through fixed defects (e.g. a pebble in a mud puddle), 
then cracks may still nucleate along the same points in each generation, but then grow to intersect with 
each other at right angles, again creating a statistically similar pattern anew in each cracking stage. If any 
memory is too strong, however, the crack pattern may repeat in each cycle without change. To evolve away 
from a rectilinear pattern new cracks must be guided, but imperfectly, by previous cracks. 

A nascent crack can be guided by a crack in an earlier cycle in a number of ways. During columnar joint 
formation each crack advance follows the crack above it, which concentrates stress in its vicinity, but can 
grow downwards in a slightly different direction [2]. In a thin layer, the previous cycle's crack may not have 
fully healed, or may have damaged its vicinity. Either case lowers the critical energy release rate G c near 
its path l2T1l . Alternatively, if the previously cracked region is more compliant (lower E\) than un-cracked 
material, then near a previous crack the strain energy release rate G will be higher, as suggested by the 
form of Eqn. [3] Finally, as with polygonal terrain, the crack may leave a depression after healing 13511621 , 
locally lowering the average height h by S. In this case the total energy spent by a unit length of crack 
extension would be reduced to G c (h - S). However, the energy released per unit length remains integrated 
throughout the region surrounding the crack of width ~ h, and will be relatively unchanged. It will thus be 
energetically favourable for a crack to be confined within the depression of the former crack. The action 
of any, or a combination, of these factors means that a crack tip advancing to maximise G - G c would be 
guided to follow any previously existing crack path. 

Changes in the order at which cracks approach and leave vertices allow for the above conditions to evolve 
a crack pattern. The case of a crack guided by a network of partially healed cracks is demonstrated in Fig. 
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[3ja-f). In this example, consider a film that cracks and heals regularly, but where the regions near previous 
cracks are more compliant than their surroundings. In the first generation, as sketched in Fig. [3Ja), one 
crack has advanced from left to right, while a later crack has curved in from below to intersect it at a right 
angle. After healing, in the next cycle of cracking, a crack is guided along the path of a previous crack. 
As it approaches the area near the previous vertex the stress field around its tip becomes asymmetric (i.e. 
the shear stress intensity factor is non-zero); as has been recently demonstrated, a crack path is attracted 
by an inclusion of more compliant material l33l . If this attraction is weak, the crack will follow close to 
its original path, and the pattern will have no opportunity to evolve. However, as shown in Fig. |3jb), the 
asymmetry may cause the crack may be attracted downwards, and curve to follow a different path to its 
first incarnation. A further crack approaching the vertex in this cycle, as in Fig. [3jc) will also be deflected 
downwards slightly, as it curves to intersect the first crack. Since variations in film stresses can only be 
felt over distances of a few film thicknesses (Eqn. [T]), these changes in crack paths will initially only be 
expected in a small region near the vertex, whose size should vary directly with /, and thus h. Within this 
region the new crack paths are deflected, and the angles at which the cracks approach the vertex will be 
more equi-angled than the original T-junction. 

If there are repeated changes in the order at which cracks approach a vertex, then the resulting changes 
in the vertex position will accumulate over many cycles, and should lead to a directed drift of the vertex in 
the direction of the crack that originally approached it at right angles. For example, if in any cycle a crack 
approaches first from the lower branch of the vertex, as sketched in Fig. [3|d), the path near the vertex will 
be unstable, in that any slight perturbation would cause it to deflect to either the left or right branch. If 
the right branch is chosen, the crack will further round the pre-existing vertex, and move it further down. 
Over many cycles, if all branches near a vertex crack first with equal frequency, the vertex will stabilise 
with three equal angles of 120°, as sketched in Fig. [3Je), by symmetry. However, as shown in Fig. [3Jf), 
the order of cracking in any particular cycle may give rise to slight deviations in the immediate vicinity of 
vertices, where cracks curve slightly to intersect each other. Alternatively, this might be noticeable by a 
damage zone, or through crack tip splittings, near the vertex. 

A pattern which evolves as it passes through space, such as columnar joints, can be thought of as the 
structure formed by stacking the evolving crack patterns described here, over many cycles, as shown in Fig. 

This model of crack ordering makes several predictions. First, that any such hexagonal pattern must 
have cracks recurring repeatedly, near older cracks, but not always opening in the same direction - and, 
equivalently, that the order at which the cracks arrive around any vertex open must change with time. 
Further, that any environment with a hexagonal crack pattern should also be able to support, and have 
previously shown, a rectilinear one. Both types of pattern should scale in the same manner, with crack 
spacings a few times the stress relaxation length / ~ h, as long as sufficient stress exists to saturate the 
spacing. In general, the hexagonal patterns must have had a chance to evolve, over many crack opening 
cycles. The reverse, that rectilinear patterns are 'young', is not necessarily true, as a crack pattern may heal 
sufficiently between openings to erase any memory, and reset the pattern to that formed in a homogeneous 
film. Finally, intermediate patterns should link the rectilinear and hexagonal forms. The effects of an 
evolving pattern should be first noticeable by small crack deflections in the vicinities of crack vertices, in a 
direction that tends to equalise the angles between cracks. Further from the vertices, in the early stages of 
evolution, cracks should still appear to be approaching as like a T-junction. 

The remaining sections of this article will show how the model conditions and predictions outlined here 
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are realised for desiccation cracks in soils, polygonal terrain, columnar joints, and eroding gypsum sand 
cements. 

IV. MUD CRACKS 

Desiccation cracking in thin elastic layers, such as clays or other granular media, is a simple model system 
that is frequently used to study crack patterns. Experiments in thin starch layers have been used to show how 
the timing information of a crack pattern is encoded in the shape of T-junctions |7 ], for example, while the 
vibrating of pastes prior to drying has been used to show how crack patterns are influenced by anisotropy in 
film properties (50l[52. In nature, as shown in Fig. [TJf)-(g), mud cracks can be found with both rectilinear 
and hexagonal cells. It has been demonstrated experimentally that mud-crack patterns evolve, when they 
are repeatedly moistened and dried l2T1l . Over only a few cycles of cracking and healing, the crack angles 
approach 120°, and verticies move. In this section, these experiments are further developed, in an attempt 
to look at how the evolution of mud-crack patterns scale with film thickness. 

Slurries of bentonite clay (Acros Organics, Bentonite K-10) were prepared by mixing one part bentonite 
with 1.9 parts deionized water (Millipore) by weight, stirring the mixture to be homogeneous, and then 
immediately pouring into large flat-bottomed glass Petri dishes (diameter 185 or 195 mm, height 15 mm). 
The slurries were dried 60 cm under a 500 W heat lamp in a well ventilated room, and were fully dry 
within 24 hours. Previous work ifTTTl shows that in this system the crack spacing is about 3 times the layer 
thickness, except in layers thinner than ~3 mm, when the critical layer thickness is approached. After each 
drying phase was complete, the dishes were imaged by a digital camera (Nikon D5100). Markings made 
around the circumference of the dishes were used as fixed points to adjust the magnification, position, and 
rotation of each image to match the configuration of the first drying's image. Typically images could be 
matched to within 2 pixels, or 0.1 mm resolution. Each dish was then rewet, by spraying with a light mist 
of deionized water, until it was visibly damp throughout the entire layer (typically a mass of water between 
1.1 and 1.3 times the mass of dry clay was added), and then set back to re-dry. 

In the first two cycles of drying and wetting, the positions of vertices were measured by hand from digital 
images, and compared in 10 dishes between 2.6 and 6.9 mm thick. Two types of vertices were investigated: 
those where cracks had deflected around a vertex, and those where cracks had formed as in the first cycle. 
Examples of both types are shown in Fig. |4ja,b). The morphology of deflected vertices was sometimes 
variable. In the thicker layers, especially, there were cases of branching or en passant ifTTl cracks near 
vertices (Fig. |4jc-d)) alongside more smoothly curving intersections (Fig. |4j e ))- F° r vertices where cracks 
were deflected off their original path in the second drying cycle, the deflection scaled linearly with the 
film thickness, as shown in Fig. |4jf). Although these displacements showed large standard deviations, the 
mean displacement over many such vertices was 14+1% of h. For the remaining vertices, there was smaller 
average deflection. However, this deflection was above the limiting resolution of 0.1 mm, and scaled with 
film thickness, at 5+1% of h. 

For 4 dishes, between 2.6 and 6.3 mm thick, the evolution of the crack pattern was monitored over 10 
wetting and drying cycles. Although it was originally intended to observe over a longer period, the spray 
bottle used to wet the clays was replaced for the II th wetting cycle, and the new spray was sufficiently 
hard to dramatically modify the pattern. For these first 10 cycles of drying and wetting, the crack patterns 
were processed automatically (Fig. [5ja-d)) by thresholding grey scale images based on local intensity and 
gradient (edge detection). These images were then skeletonized, and vertices found as the branch-points 
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FIG. 4: Mud-crack patterns change between (a) their first and (b) second dryings. In (b), two vertices have changed 
the order in which cracks have approached a vertex, and shifted their position significantly as a result (their positions in 
the first cycle are shown by blue points), while two vertices (red points) have remained relatively undistorted. In thicker 
layers, branching of cracks is sometimes noticed, for example between (c) the first and (d) second drying, alongside 
(e) smoother vertex deflections, (f) The average displacement of vertices with cracks that have curved away from their 
original position (blue discs) and vertices which have not (red circles), are both consistent with a linear dependence on 
layer thickness. In (a)-(e) scale bars are 1 mm. 

of the thresholded image skeleton. Using a Matlab-based particle tracking code lfT2l . the motion of the 
vertices was followed. As with previous experiments [21 ], cracks could sometimes disappear for several 
cycles, only to reappear close to their original location. For this analysis, only vertices that were present in 
all 10 crack patterns were analysed. After this time the average vertex has moved between about 0.5 to 1 
mm, depending on the layer thickness. As shown in Fig. |5Je), when scaled with the layer thickness h, the 
motion from all patterns collapse onto a single curve, which also overlaps data published previously ETI 
for bentonite clay 2.9 mm thick. In all cases, there was a significant change in the average vertex positions 
within the first 3 or 4 cycles, followed by a much more gradual vertex motion. Previously l2TTl . it was 
observed that 25-35% of vertices changed order of formation in any particular generation. The slowdown 
in vertex motion thus appears to be identifiable with the time it takes for all vertices to change order of 
opening once. Furthermore, the average displacement after 4 cycles, 15+1% of h, agrees with the vertex 
shift, shown in Fig. |4jf), for the subset of cracks which changed the order of arrival at a vertex between the 
first and second drying cycles. 

The angle at which cracks approached vertices was also measured, automatically. A circle of small, but 
arbitrary, radius was drawn around each vertex (branch-point) of the skeletonized image, the points at which 
the skeleton crossed the circle were found, and the angle between crossings noted. A radius of about 30 
pixels (2 mm) was typically a good compromise between the errors associated with too small a sampling 
radius, and the desire to measure the immediate vicinity of the crack vertices; for small radii, a systematic 
underestimation of the largest (near 180°) angle between cracks was noticed. As shown in Fig. |6j there is a 
gradual but systematic equalization of all crack angles towards 120° over many cracking cycles. 

Finally, Fig. [5jf) shows that the direction of the displacement of vertices between the first and tenth 
drying cycles is well-correlated with the direction of the crack that originally approached the junction at 
the smallest angles to its neighbours. In other words, the vertex is moving preferentially along the direction 
required to stretch a T-junction into a Y-junction, the direction of the crack that originally approached the 
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FIG. 5: The motion of vertices can be measured over many drying cycles. During analysis (a) a raw image is (b) 
thresholded, (c) skeletonized, and then (d) vertex positions are found as the branch-points of the skeleton. Shown here 
are the positions of three vertices over ten drying cycles, graded in colour from the first (blue) to the tenth (red) cycle, 
(e) The average displacement of all vertices increases with the number of drying cycles, and also scales with the layer 
thickness. Representative error bars are given for h = 3.6 mm. The data from a previous study ['21 1 where h = 2.9 
mm are also shown. The changes are most apparent in the first four cycles, and then slows down to a gradual drift. 
The arrow shows the measured vertex displacement when cracks round a vertex, from Fig. [4jf), f° r comparison. The 
direction (f) in which the vertex travels is well-correlated with the direction of the crack that originally approached 
the vertex with the smallest joint angles. The drift angle 6 is the difference between the vertex displacement measured 
between the first and tenth cracking cycles, and the angle of approach of this crack, as shown in the inset. Data here are 
shown for the 2.6 mm thick sample. 
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FIG. 6: The angles at which cracks approach vertices in drying mud-crack patterns move toward the 120° angles of 
Y-junctions over a few cycles of cracking and healing. Shown are the probability distribution functions for the angles 
between cracks measured at radii of 1.7 (red), 2.0 (black - with representative error bars) and 2.7 (green) mm from the 
vertices of a 2.6 mm thick dried bentonite layer, after its (a) first, (b) fifth, and (c) tenth drying is complete. 

vertex at near two right angles. 

V. POLYGONAL TERRAIN 

Polygonal terrain consists of networks of thermal contraction cracks, typically a few tens of meters apart, 
which form in ice-cemented soils (35). It is common in terrestrial circumpolar regions. The thermal stresses 
arising from the intense winters of these regions are sufficient to open cracks a few cm wide, and several 
meters deep (35j[38]|4l]]. While open loose soil, snow, and wind-born detritus may fall into the crack. Thus 
the crack does not heal completely when thawed, and is likely to reopen in the following winter. These 
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FIG. 7: Polygonal terrain in the Dry Valleys, Antarctica, (a) Oblique/arial view of cracks near New Harbor delta, 
of age 1 - 2 ka |62|. Crack positions are highlighted by white arrows along the perimeter of the image. The vertices 
shown are distorted, and appear to be deflected in the direction indicated by black arrows from what would otherwise 
be T-junctions, (b) Close-up image of a vertex of a Y-junction that shows a small 'core' where cracks in any particular 
year can interact to approach each other at right angles, (c) Black and Berg installed pairs of rods here bracketing the 
crack pattern in 1963. By 2007 this crack (highlighted by white arrows) has shifted to lie outside the rod pair, (d) At 
Victoria Valley, where the terrain is 10-12 ka |62|, the pattern is more hexagonal. 

patterns have analogies to the evolving mud crack patterns, where each seasonal thermal cycle is equivalent 
to one cycle of drying and wetting. However, the net addition of material into the crack creates a wedge 
of ice or sand that grows at rates of order a millimetre per year and can lead to a slow turnover of the 
near-surface soil (4l l42ll62) . 

Recently, perhaps as a result of the landings of Phoenix [44 ] and Curiosity (26) on Martian polygonal 
terrain, much work has gone into analysing polygonal crack patterns via satellite imaging. Large datasets 
have been collected for sites on Earth and Mars (221 ESI E3 E2). These databases have been used to 
demonstrate that polygonal terrain shows statistical trends such as are also found in foam networks (57) and 
that the vertex positions contain spatial correlations [ 27 ], and work characterising these networks is ongoing 
and interesting. Comparisons with ground observations have shown, however, that very high-resolution 
images (~0.25m/pixel) are required to correctly measure the features of polygon vertices (271 [56), making 
detailed interpretation of satellite images difficult. 

More precise results have been gathered through a small number of very long-term field experiments 
into the dynamics of polygonal terrain. In Arctic Canada, Mackay installed cables across the positions of 
ice-wedges, and measured the time of crack opening as the time at which the cables broke, during several 
consecutive years [38-42]. At Illisarvik a small permafrost lake was drained, and the development of new 
ice- wedges was studied over the subsequent twenty years (39l[4T][42l. The initial crack pattern was strongly 
influenced by local features, but the most common crack junction was a rectilinear T-junction (42) . There 
were rapid changes in the growth of ice wedges within the first few years, before ceasing due to vegetation 
coverage (42). On Garry Island, several well-established crack patterns were instrumented (38). 

Not all cracks opened each year, yet cracks could reappear in the same spot after several-year absences. 
This is similar to drying experiments in clays, where cracks sometimes reappeared after a hiatus of several 
drying cycles. The cracks themselves appear to initiate between the position of the buried ice- wedge, and 
the soil surface (39] 02) . The interpretation of this is that the wedge itself acts as a memory in this system, 
as it is preserved over many generations, and is a foreign material that is weaker, or more compliant, than 
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the surrounding soil (35l[39l[42). The cracks advanced slowly, in many cases taking several days to travel 
from one vertex to the next [40 1. It was noted that the cracks could shift positions slightly, by several cm, 
in different years (38). The cracks often showed en passant style doubled paths 1401 , similar to the example 
mud cracks shown in Fig. [4jd), and other branching structures. It was also seen that cracks initiated at 
different points, and opened in differing orders, in different years [40 1. 

Another long-term experiment was initiated by Black and Berg in the Dry Valleys of Antarctica (U[6), 
and developed further by Sletten et. ah [62]. At several sites, pairs of metal rods were driven into the 
permafrost, bracketing the cracks of polygonal terrain, and the separation of these rods carefully measured. 
An intermittent series of measurements has been made on rod separation spanning 40 years (U |6] [62), 
showing that wedges were growing at rates of -0.5 mm/year. The sites developed include a range of ages at 
New Harbor/Taylor Valley (1-2 ka, Figs. [TJd) and^a-c)), Victoria Valley (10-12 ka, Fig. [TJd)), and Beacon 
Valley (>1 Ma, Fig. [TJe)), which also display different patterns (62). In older patterns the crack pattern 
is dominated by Y-junctions, while at New Harbour it is closer to a rectilinear patten [62]. Based on the 
relative ages of these sites, and their degree of ordering, it has been suggested that an initially rectilinear 
pattern of polygonal terrain cracks evolves towards a more hexagonal pattern over several thousands of 
years [62]. 

In the vicinity of the New Harbour site, cracks show additional features that agree with the means of 
evolution proposed here. As shown in Fig. [7ja), there are regions around many vertices where the cracks 
curve suddenly to meet the vertex positions, as was sketched in Fig. |2jc,d), or shown in mud cracks in Fig. 
|4jb,e). This suggests that the vertices have moved in the direction perpendicular to the cracks that begin 
their approach at -180°, away from their original positions. In this case the motion should be recorded in 
an observable asymmetry of the wedge material beneath the crack pattern; the deposits on the side that was 
originally the primary branches of a T-junction should be thicker than the deposits on the side that originally 
had two right angles. This prediction could be checked by excavating a statistically significant number of 
vertices. 

Some vertices showed a complex inner structure, that appears to record how active crack branches curving 
to intersect each other slightly with different orders of opening in different years, as shown in Fig. [7jb). 
Finally, in December 2007, at Black and Berg's New Harbour site, it was noted that cracks at 2 of the 44 
pairs of rods installed in 1962 had shifted position to now lie outside the pairs of marker rods (see Fig. [TJc)). 
As with Mackay's more direct measurements of crack positions in Arctic polygons, this indicates that the 
crack positions can change over time. 

Taken together, the observations of polygonal terrain show that the fracture network is persistent, but that 
the pattern does change over time. The persistence is related to the fact that cracks leave permanent records 
of their existence within the soil, which preferentially nucleate and guide later generations of cracks. The 
order of crack opening must also vary, as not all cracks are active in each year. All these features were seen 
evolving mud crack patterns, and were essential in developing a model of crack pattern evolution. 

VI. COLUMNAR JOINTS 

Columnar joints in lava are a well-known geophysical pattern, and have attracted scientific curiosity for 
centuries [8,31]. The same pattern can be observed in other media, such as quenched glass (18), thermally- 
shocked sandstone (60), and vitrified ice (45) . It also appears in desiccated starch, where it has been 
partiularly well-studied OfT^l22l^inilM^l5l[6^ . The columns in starch are statistically identical 
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FIG. 8: Columnar joints in lava, (a) A cross-section of a colonnade near Fingal's Cave, Staffa, Scotland appears 
well-ordered, but (b-d) cracks very near the Y-junctions curve slightly, reflecting interactions during a particular growth 
increment, (e) A single column near Whistler, Canada, showing lateral striae bands. These can be used to infer the 
crack front growth velocity, which is faster (the striae are smaller) near the base of the flow, (f) The crack front velocity 
in basaltic (black) and other lavas (green) is inversely correlated with the size of the columns, (g) as it is also in starch 
for constant (green) and slowly varying (black) front velocities (panels (f,g) adapted from (23)). Curves show the 
expected correlation for constant Peclet numbers, (h) Plumose, or hairlike marks on striae, here sketched over with 
chalk, indicate the direction of crack propagation of each crack increment, (i) The crack growth direction recorded on 
striae can be toward (white) or away (grey) from any vertex. Shown here are the direction of crack propagation for 
sequence of striae on the faces of three columns, near different vertices. 



to columns in lava (241125) . They are, however, usually 0.1-1 cm in diameter, while columns in lava are 
typically of order 10 - 100 cm across. 

In lava, columnar joints form by the gradual extension of thermal contraction cracks into a cooling lava 
body, with the cracks growing perpendicular to the isothermal surfaces at the time of cracking (43) . Thus, 
at any particular time, the crack tips represent a network in a thin region near the surface defined by the 
glass transition temperature isotherm (59), and the columnar forms can be interpreted as the outlines of 
this network as it evolved in time, whilst propagating through space. Columns may grow in any direction, 
and in many cases form both an upper and lower colonnade in a single flow unit, as the lava cools from 
both upper and lower boundaries (eg. O [37)). In starch, columns grow normal to surfaces of constant 
moisture content (20|[23), and track the position of the so-called funicular-pendular transition, where the 
water transport mechanisms of capillary flow and vapour diffusion cross over in efficiency (20). 

The crack pattern at the surface of lavas observed on Hawaii are rectilinear in appearance, and dominated 
by T-junctions El [55). The evolution of column vertices was observed at the Boiling Pots, Hawaii, where 
columns with an average size of ~ 20 cm 11131 were seen to develop from T-junctions to Y-junctions, over a 
distance of 1 m (2). Under careful inspection most vertices were classified there as 'pseudo-Y' junctions, as 
the cracks curved slightly in the immediate vicinity of the vertex [ 2 ], similar to the sketch here in Fig. J3jf). 
This small deflection was used to argue that all cracks did not interact at a vertex synchronously (2). Similar 
vertices can be observed in many colonnades; an example from Staffa is shown in Fig. [SJa-d). In starch 



13 



colonnades, an evolution from a disordered surface pattern to a hexagonal mature pattern is also observed 
l24l . For columns 1-2 mm across, the pattern orders over a distance of ~1 cm from the drying surface l24l : 
this ordering region can be seen in Fig. [TJb). Although these observations are sparse, it seems that in both 
starch and lava the hexagonal columnar pattern is well ordered after 5-10 column diameters of growth. 

As cracks advance in lava, they do not do so continuously, but intermittently, leaving striations on the 
sides of columns as shown in Fig. [8je,h). Each stria represents a single growth increment: once a crack tip 
activates it rapidly advances through brittle lava, following laterally near the path of the previous crack's tip. 
In the forward direction of column growth, however, the tip blunts and stops as it intrudes into a warmer, 
more plastic zone, and each crack advance thus leaves a broad thin mark on the side of a column [58]. 
These striae record the cooling conditions of the lava at the time of their formation, and can be interpreted 
accordingly (2[|T4l[25l|5U. In particular, the laterally-averaged height of the striae represents the distance 
between two isotherms in the cooling lava body, and can be used to measure the cooling rate at the time of 
fracture l25l . 

During the cooling of a lava flow, water or moisture in the environment commonly acts an efficient 
transporter of heat, and the thermal transport problem can be approximated as that of an arbitrarily thick 
layer, cooled entirely to 100°C, bound to a purely diffusive region over which the temperature transits to that 
of the internal melt (l5][25][28]l. The majority of the thermal gradient is confined to this thin layer, with a 
thermal diffusivity D, and all isotherms traveling together at some velocity v. This linear advection-diffusion 
problem has a steady- state temperature profile T, that follows 

V 2 r + Pe^=0 (5) 

oz 

where Pe = vL/D is the Peclet number, and L is a characteristic length proportional to the thickness of 
the diffusive layer. The velocity v can be inferred, quantitatively, from the average striae heights along any 
particular cross-section of a colonnade l25l . As the thermal gradients are confined to the advective-diffusive 
layer, so are the elastic stresses, and according to the arguments presented in Section 2, the crack spacing 
should be proportional to this layer thickness, and inversely proportional to v [25|. As shown in Fig. [8jf), 
this is the case: large columns cooled slower than small ones. 

In starch, the moisture concentration is functionally equivalent to the temperature, in generating stress. 
There is a bottleneck in water transport, however, as the two mechanisms of vapour diffusion and capillary 
flow cross over in efficiency at some critical moisture concentration [20]. Most of the stress gradients in a 
drying starch body are confined into a thin layer, near this bottleneck, where an advection-diffusion problem 
can be posed for the moisture concentration, 

V 2 + Pe h ^ =0. (6) 

oz 

Here D h is the minimum in the moisture diffusivity and Pe/j = vL/Dh E3 E3- As shown in Fig. [SJg), the 
size of columnar joints in starch is inversely proportional to their cracking front velocity v, and when scaled 
by D h are found at the same Peclet number as those in lava [ 23 ]. 

Columns can evolve by slight variations of the crack plane in each crack increment fT4lL or by individual 
cracks overshooting or rounding a vertex formed in the previous cracking cracking cycle (2). In general, 
each crack increment is guided by its predecessor, which concentrates stress along its tip, but can advance in 
a slightly different direction. In starch this evolution has been followed through scanning x-ray tomography 
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fl3\ l22l l24l |48l and magnetic resonance imaging (46). The topological changes of the pattern can be 
described by three cases: two vertices may collide and exchange position; crack faces can fail to propagate, 
merging two (or more) columns; and new columns can initiate from existing vertices [22]. For neighbouring 
columns of different sizes, there is also a strong tendency to equalise the column area l22l . suggesting long- 
term dynamics that are not covered by the model presented here. 

Finally, cracks in columnar lava do not open in the same order, in different cracking increments. The 
striae markings record the direction of crack propagation of each crack increment in their so-called plumose 
structure fT4l . From this, one can tell whether a crack advanced toward, or away from, any vertex in any 
particular growth increment. As has been seen elsewhere (21 [141 [T9lL and as shown in Fig. [SJh), this 
direction can change along subsequent striae. At one colonnade near Bridge Lake, a site described in Ref. 
|[T9lL the directions of growth of 835 striae were measured. Some example growth direction sequences are 
shown in Fig. [8ji). It was found that 64+2% of the time, one striae followed the same direction as it's 
predecessor. Thus, while each striae direction is not randomly chosen, the sequence at which cracks are 
activated at any vertex changes frequently during a column's formation. 

Columnar joints develop in a very similar way to the other hexagonal crack patterns described and mod- 
elled here, where the cycles of cracking and healing have become a stacked set of crack advances. In each 
advance cracks are guided by the cracks that preceed them, but are allowed to take a slightly different path. 
The sequence of cracking at any vertex is variable, however, and the vertices themselves interact and move 
over time. The crack pattern rapidly evolves from a rectilinear, T-junction dominated pattern to a hexagonal 
one, but retains a slight asymmetry in the immediate vicinity of vertices that is indicative of the order of 
crack opening in any particular growth increment. 

VII. ERODING CRACK PATTERNS 

Recently, a crack pattern that forms on cemented gypsum sands, on the windward (and therefore eroding) 
faces of dunes, has been described OHO), and will be briefly reviewed here. These crack patterns, shown 
in Fig. [IJc), are found on the dunes of the White Sands National Monument, in New Mexico. Genetically, 
they have been argued to be similar to columnar joints, but where the cracks are advancing into the dunes 
as the cemented sand surface is removed ifTOl . Crack depths of 0.1 cm to 45 cm were reported, and polygon 
sizes from 5 to greater than 70 cm. As with naturally-dried desiccation cracks in starch (24l[48) and thermal 
contraction cracks in coke [34], the polygon size was noted to increase with depth. 

It was found that cracks changed directions abruptly as vertices were approached ifTOli . In this, they 
appear similar to the mud crack experiments reported here and in Ref. [21], and to the devloping polygonal 
terrain at New Harbour (62). The angles between cracks at the junctions were measured at 4 sites, and 
found to deviate significantly from those of orthogonal T-junctions IfTOl . Only in very thin (1 mm) cements, 
or cases where the cracking directions were strongly influenced by the dune slope, were T-junctions typical. 

A convincing model for the origin of these cracks was also proposed 1 10]. The cycling of moisture in 
the dunes allows for gypsum to dissolve and re-precipitate near the dune surfaces, cementing the sand. The 
same water, as it evaporates, causes stresses which fragment the cemented layer. Finally, wind erodes the 
stoss face of the dune, allowing the wetting (and cracking) layer to slowly intrude into the dune. 

Several features of these patterns, such as the local changes in crack orientation near vertices, and the 
opportunity for dynamics as the face erodes, are consistent with an evolution similar to that outlined here. 
However, direct evidence for this is limited. At one site 13 junctions were excavated, and measured at the 
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surface, 2 cm, and 10 cm depth. Although it was stated that cracks became straighter, and more ordered, 
with depth, the distribution of angles from this sample did not significantly change ifTOli . As it is an inter- 
esting pattern, potentially related to ordered cracks and water cycling on Mars |9], additional work is to be 
encouraged. 

VIII. SUMMARY 

The crack patterns discussed here, such as are commonly seen in dried mud, frozen soils, or old lava 
flows, can take the appearance of a rectilinear network, with cracks intersecting preferentially at near right- 
angled T-junctions, or a hexagonal network of equi-angled Y-j unctions, or sometimes an intermediate state. 
A genetic model was sketched which could account for these observations, amongst others, as the result of 
simple interactions. A crack curves to maximise the difference between the strain energy release rate, at the 
moment of cracking, and the cost of creating the new crack surfaces - an established principle of fracture 
mechanics l36ll . This is a local energy consideration: the crack tip only senses its immediate environment, 
and responds accordingly. When a layer repeatedly cracks and heals, the positions of the previous cracks 
may be weaker, or more compliant, than the surrounding material. Similarly, in a crack network like that of 
columnar joints, that is growing through space, the stress state in an un-cracked layer is strongly influenced 
by the prior pattern. In this way cracks are guided along the paths taken in previous cracking cycles. At 
vertices, however, this argument introduces asymmetries, and a growing crack can be deflected. If the order 
in which cracks arrive at a vertex changes with time, in different cracking cycles, these deflections can 
accumulate and will distort the vertex, moving it in a direction to equalise the angles between all cracks, 
creating a Y-junction. The requirements and predictions of this process were outlined in Sections 2 and 3. 

This model was compared to and elaborated by presentations of evolving crack patterns in desiccated 
bentonite clays, polygonal terrain in permafrost, columnar joints in lava and starch, and a pattern of desic- 
cation cracks that evolves as it erodes. These patterns were found to share similar features, such as a scaling 
of the crack spacing with the layer thickness, changes in the order of crack opening, and distortions of the 
crack pattern in the immediate vicinity of vertices. Furthermore, it was shown in bentonite clay layers that 
were repeatedly dried and rewetted, that the drift of vertices away from their original positions, during this 
evolution, scales with the thickness of the cracking layer. This rate may may also depend on the materials 
involved - as permafrost apparently takes much longer to order than a few tens of cracking cycles, mea- 
sured there by years. The direction of this motion was seen to be precisely that predicted to deform vertices 
towards a hexagonal crack pattern. 
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